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The role of magnetic and electric perturbations in the quantum Lifshitz model in 2 + 1 dimensions 
are examined in this paper. The quantum Lifshitz model is an effective field theory for quantum 
multicritical systems, that include generalized 2D quantum dimer models in bipartite lattices and 
their generalizations. It describes a class of quantum phase transitions between ordered and topolog- 
ical phases in 2 + 1 dimensions. Magnetic perturbations break the dimer conservation law. Electric 
excitations, whose condensation lead to ordered phases, have been studied extensively both in the 
classical 3D model and in the quantum 2D model. The role of magnetic vortex excitations whose 
condensation drive these systems into a Z2 topological phase has been largely ignored. Recent nu- 
merical studies claim that the quantum theory has a peculiar feature: the dynamical exponent z 
flows continuously and the quantum theory is hence unstable to magnetic vortices. To study the 
interplay of both excitations, we perform a perturbative renormalization group study to one loop 
order and study the stability of the theory away from quantum multicriticality. This is done by 
generalizing the operator product expansion to anisotropic models. It is found that the dynamical 
exponent does not appear to flow, in contrast to the classical Monte Carlo study. Possible reasons 
for this difference are discussed at length. 

PACS numbers: 64.60.Ht, 71.0.Hf 



I. INTRODUCTION 

In this paper we examine the scaling behavior of the 
dynamics of quantum dimer models near quantum criti- 
cality and their generalizations. Quantum dimer models 
are 2D strongly correlated systems with a rich phase dia- 
gram that includes both ordered and topological phases. 
The quantum dimer model was introduced by Rokhsar 
and Kivelson (RK) in 198& 1 as a model to describe spin 
liquid states within a large spin gap. Such quantum dis- 
ordered phases presumably arise in sufficiently frustrated 
quantum antiferromagnets so that all spin-ordered states 
are suppressedP^ Typical examples of frustrated sys- 
tems that can be described by quantum dimer models 
include antiferromagnets in transition metal oxides on 
Kagome and pyrochlore lattices.®^ 

The most interesting feature of the quantum dimer 
models is the rich variety of phases that are accessible 
from a quantum multicritical point with special proper- 
ties. The phase diagrams of these models are known by 
a combination of strong coupling arguments and from 
the existence of a special choice of parameters, know as 
the Rokhsar-Kivelson (RK) point, at which the quantum 
Hamiltonians take the form of a sum of local projection 
operators which allows for the determination of their ex- 
act ground state wave function (at the RK point only). 1 
The resulting ground state wave functions have the form 
of a liner superposition of dimer configurations with lo- 
cal weights. A consequence of this structure is that the 
norm of the ground state wave functions of these models 
(at their RK points) is equal to the partition function of 
classical dimer models on the same 2D lattice whose (lo- 



cal) Gibbs weights are the square of the amplitudes of the 
wave function. For the simplest quantum dimer model, 
the ground state wave function is equ al to the short-range 
resonating valence bond statepSEI] Because of the close 
connection between the ground state wave function of 
generalized quantum dimer models and two dimensional 
classical statistical mechanical systems, a great deal of 
information about the static properties of such models 
are known. Given the knowledge of the exact ground 
state wave function, it is also possible to compute the 
equal-time-correlation functions by mapping them to a 
2D classical counterpart. It has also been instrumental 
in exact computations of the entanglement entropy in 
two dimensions j 12 " 16 ! thus extending the known results 
for one-dimensional systems.^^Sl 

An interesting system of this type is the quantum 
eight-vertex model of Ref.[3T]. This model is a gener- 
alization of the 2D quantum dimer model. In this case 
the degrees of freedom are arrows placed on bonds of the 
square lattice. The allowed configurations are shown in 
Fig. [T] This model also has an RK-type point at which 
the ground state wave function has local weights with 
the same form as the Gibbs weights in the corresponding 
Baxter (or eight- vertex) model on the square lattice, i.e. 
the Baxter wave function. The partition function for the 
classical model (which in the present context is the norm 
of the Baxter wave function) was solved by Baxter P^l The 
classical Baxter model has ordered and disordered phases 
separated by two critical lines (with continuously varying 
critical exponents) that meet at a multicritical point (see 
Figure [2} P 1 

Ardonne et alW^ used the known results of the 2D clas- 
sical Baxter model to determine the phase diagram of the 
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FIG. 1. The amplitude of the wave function depends on the 
type of vertices that enter in each configuration. Vertices 
in the (a) six- vertex model with amplitudes ("fugacities") a, 
b, and c, respectively, and (b) the additional vertex in the 
eight-vertex model with fugacity d. These vertices and their 
inverses with all arrows flipped, form the six vertex and eight 
vertex models respectively. The heights are displayed in quad- 
rant. To do so, one arbitrarily picks a starting plaquette and 
moving in a counter clockwise direction, raises the height vari- 
able by 1 if an arrow going into the vertex is crossed or de- 
creases the height variable by one if an arrow going out of the 
vertex is crossed. The configurations in |(b)| break the ice rule 
of the six-vertex model, and are vortices of the height field. 

quantum counterpart at its RK "point" , shown in Fig. [2] 
as well as the behavior of the equal-time correlators in 
the phases and at (quantum) criticality. A key result 
of that work is that the disordered phase of the classi- 
cal model is the Z 2 topological phase of the quantum 
model. In particular there is a special point deep inside 
this phase (where all the weights are equal) at which the 
quantum Hamiltonian is identical to the Hamiltonian of 
Kitaev's toric codef^ Hence, the Baxter wave function 
can be used to study a quantum phase transition from 
an ordered phase to a Z 2 topological phase by crossing 
one of the critical lines. 

However, the knowledge of the ground state wave func- 
tion (and hence of all the equal-time correlation func- 
tions) is not sufficient to determine the universality class 
of the quantum phase transition. We therefore have the 
strange situation where much is known about the phases 
on either side of a critical point, but the theory at the 
critical point is incomplete. Although the quantum dimer 
models (and the quantum eight-vertex model) are not 
exactly solvable, variational arguments by Rokhsar and 
KivelsorP suggested that the dynamic critical exponent 
of these problems is z = 2. This observation was for- 
malized by Henley^ and in more detail by Moessner et 
aZpO who suggested an effective field theory, obtained by 
coarse-graining the (dual) height representation, which 
has z = 2. Ardonne et alW^ studied the consequences 
of this theory, which they dubbed the "quantum Lifshitz 
model." They conjectured that this model represents the 
universality class of generalized quantum dimer models 
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FIG. 2. Phase diagram of the quantum eight vertex model 
(from Ref. [21]). Bold lines represent quantum critical lines. 
The effective field theory along the lines labelled "six-vertex 
model" is the quantum Lifshitz model. Here c and d represent 
amplitudes of the wave function (see text). The quantum dis- 
ordered phase is a Z2 topological phase and the point labelled 
"Kitaev" is the toric code. At the level of the wave function, 
the two critical lines are related by Kramers- Wannier duality. 
We will be studying the theory near the multicritical point, 
X , into the full eight vertex model along the critical line with 
d 2 ,c 2 £ 0. 



and both critical lines of the quantum eight-vertex model. 

The main purpose of this paper is to extend the anal- 
ysis of Refs.[57] and [55] to include the effects of such 
magnetic (vortex) perturbations on the quantum Lifshitz 
model and to assess their effects on the quantum uni- 
versality class and more specifically on the dynamical 
scaling exponent z. We will specifically look at the case 
of the quantum eight-vertex model. Although the Bax- 
ter wave function obeys KW self-duality, the quantum 
Hamiltonian of the quantum-eight vertex model (whose 
exact ground state is the Baxter wave function) does not. 
This implies knowing the static properties of the corre- 
lators along the six-vertex model line it is possible to 
determine what they are along the dual line, i.e. the line 
d 2 ,c 2 ^ in Figure [2j using KW duality. However, since 
the full quantum Hamiltonian is not invariant under KW 
duality, we will not be able to infer the properties of the 
dynamical correlators along the dual line from their be- 
havior along the six-vertex model line. In particular it is 
possible that the dynamical exponent z may not be equal 
to 2 along the dual line. More generally, we will show in 
this paper that along the dual line, the effective field 
theory (the quantum Lifshitz model) is in general unsta- 
ble under the action of perturbations that break the RK 
condition. Thus the dynamics of this system with a RK 
condition is very different than without it. 

Our strategy is to consider the quantum stability of 
the quantum Lifshitz model under the effects of generic 
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local perturbations and to use the renormalization group 
to infer the behavior of the system. To this end we 
include in the effective Lagrangian in 2+1 dimensions 
operators which describe the local creation and destruc- 
tion (in space and time) of electric and magnetic charges, 
i.e. local fluctuations of dimers (or arrows in the Bax- 
ter language) and processes involving "holons" (or equiv- 
alently that break the symmetry of the Baxter model 
down to Z2). These local perturbations violate the RK 
condition. The recent Monte Carlo results of Isakov 
and coworkers,^ which motivated our work, simulate 
the dynamics of a system in which the RK condition is 
obeyed. Marginally irrelevant perturbations and relevant 
operators which move the system away from quantum 
criticality are also included in the effective Lagrangian. 
The RG is constructed using an operator product ex- 
pansion (OPE) adapted to anisotropic systems (with dy- 
namical exponent z > 1.) Consistently with the earlier 



rcsult a * ' (and with the Monte Carlo simulations 2 ^!) we 
find that the z — 2 dynamics of the quantum dimer model 
(and the quantum six-vertex model) is (marginally) sta- 
ble in the absence of "magnetic" (vortex) perturbations, 
but it becomes unstable if these perturbations are in- 
cluded. In this case the actual behavior of the system 
cannot be determined perturbatively. We conclude that 
the likely result is that either the (quantum) transition 
becomes first order or, if it is continuous, it is con- 
trolled by a non-perturbative fixed point (presumably 
with "generic" z = 1 dynamics.) 

This paper is organized as follows. In Section [XT] we 
summarize the current status of the problem of the quan- 
tum dynamics of generalized quantum dimer models in 
2D. In Section |III| we describe the quantum Lifshitz 
model in detail. In Section |IV| we present the one loop 
RG calculation performed, details of which (including 
the OPE) can be found in Appendix B] by the interested 
reader. In Section |V] we discuss the differences which re- 
sult from dynamics that preserve the RK condition (and 
in the Monte Carlo simulations of Rcf. [29 and more gen- 
eral dynamics resulting by perturbing the quantum Lif- 
shitz model as we use here. Finally, in Section |VI| we 
discuss our results and some open problems. 



II. QUANTUM DYNAMICS OF GENERALIZED 
QUANTUM DIMER MODELS 

The validity of the quantum Lifshitz model of the ef- 
fective field theory of generalized quantum dimer mod- 
els at their RK points was established in Refs. [3T] and 
126] . Whether this theory describes the quantum dynam- 
ics away from RK points is far less obvious. In this sec- 
tion we summarize what is known about this issue. 

The validity of this conjecture has been considered 
by several authors. The combined effects of perturba- 
tions due to so-called "charge" operators (defined in later 
sections) and various non-linear gradients of the coarse- 
grained height field this quantum phase transition was 



studied using perturbative renormalization group (RG) 
methods! 27 * 28 ^ It was found that in many cases the RG 
flows destabilize the quantum Lifshitz fixed line, whether 
by rendering the quantum phase transition first order or 
by replacing it with an infinite sequence of "tilting" tran- 
sitions (dubbed an incomplete devils' staircase). Still in 
some cases the the quantum Lifshitz fixed line was found 
to be perturbatively stable (up to marginally irrelevant 
operators). This is what happens in the quantum dimer 
model on the square lattice^ and in the quantum eight- 
vertex model a long t he six-vertex model line (the d = 
axis of Fig. [2}) 21 | 31| E^ Nevertheless, these perturbative 
RG studies cannot determine (or exclude) the possibility 
that one unstable trajectory of the RG flow may drive 
the system into another fixed point (not perturbatively 
accessible) which would be, presumably, in the conven- 
tional z = 1 universality class. 

It is known that excitations that break the dimer con- 
straints (by turn a bipartite lattice into a non-bipartite 
lattice)^] or that break the ice-rule constraint of the six- 
vertex model (by addition of vertices with weight d in 
the eight- vertex models (see Figjl])) trigger a transition 
into a Z2 topological phase. In terms of the dual height 
models (and hence of the quantum Lifshitz model) these 
perturbations are vortex (and hence "magnetic") opera- 
tors. The effects of these magnetic perturbations on the 
quantum criticality of these systems is presently poorly 
understood. What we know is that in the case of the 
quantum dimer models on the square lattice perturba- 
tions that break the bipartite property of the lattice are 
relevant a nd dr ive the system immediately into a topolog- 
ical phaseP32l The same behavior is found along the six- 
vertex model line of the Baxter wave function!^! Using 
the Kramers- Wannier (KW) self-duality that the Bax- 
ter wave function has, Ardonnc et al. were also able 
to compute the scaling behavior of the excitations along 
the "dual" critical line, the phase boundary between the 
topological phase and the ordered phase, along the criti- 
cal line where d 2 7^ and c 2 7^ (see Figure [2J. However 
these results do not determine what effects, if any, these 
excitations have on the dynamics and its scaling behav- 
ior. 

While perturbations of the quantum Lifshitz theory 
has been a topic of interest to many, these studies so far 
have been limited to "charged" excitations. Physically 
these perturbations can be thought of as point defects 
in the height model and are perturbations that preserve 
the ice rule (two arrows in and two arrows out at each 
lattice site) of the six vertex model (see Figure [T]) . As 
such they are limited to the six-vertex limit of the full 
quantum eight vertex model. In the full quantum theory, 
magnetic vortex like excitations are also allowed. These 
correspond to topological defects to the height represen- 
tation of the model. On the lattice, they are violations of 
the ice rule bayed by the allowed configurations (shown 
Figure [T]). Whereas "charge" excitations can be handled 
by standard techniques, the effect of magnetic vortex ex- 
citations is less trivial and because of this, much is not 
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known in the full quantum eight vertex model. It is the 
effect of these magnetic vortex excitations with the pres- 
ence of charged excitations that forms the subject of this 
paper. We will examine the effect of these perturbations 
along the critical line within the full eight vertex model 
(see Figure [2]). 

Recent numerical studies by Isakov et. alP^have given 
further impetus for the study of the effects of magnetic 
excitations in the quantum eight-vertex model. In their 
study, they used a stochastic classical approach to simu- 
late the classical critical dynamics of the classical eight- 
vertex model. In this approach one does not directly 
simulate the dynamics of a quantum system in 2 + 1 di- 
mensions but instead, relies on strictly enforcing the RK 
condition, and does a classical simulation of the relax- 
ation to equilibrium of a 2D system. In this approach 
the effects of magnetic excitations is taken into account 
while requiring the system to strictly obey the RK con- 
straint. Indeed, the classical stochastic theory is essen- 
tially a Langevin equation (or equivalently a master equa- 
tion) for the classical stochastic time evolution, and has 
a path integral representation developed by Martin, Sig- 
gia and RoseP^ This theory is known to have a form 
of supersymmetrySsHSU which ensures that the Langevin 
dynamics is strictly obeyed. In contrast the quantum Lif- 
shitz model is defined by a quantum mechanical action 
in 2 + 1 dimensions. 

At the level of the free quantum Lifshitz models (i.e. 
along the six vertex line) the stochastic classical simula- 
tion and the quantum mechanical path integral are essen- 
tially equivalent to each otherPS However, they do not 
necessarily agree once the perturbations are included. In- 
deed in the Langevin approach one perturbs the wave 
function, whereas in the 2 + 1 dimensional quantum 
field theory one perturbs the full quantum Hamiltonian 
whose form is no longer the sum of projection opera- 
tors and, hence, in general is not equivalent to a classical 
Langevin dynamics problem. Isakov and coworkers car- 
ried out a detailed classical Langevin simulation of the 
eight-vertex model wave function. Along the six-vertex 
model line, d = 0, their results are consistent with z = 2 
along the e ntire fixe d line, in agreement with earlier ana- 
lytic results^* 28 * 30 ! and with the results that we report 
here. However, for d > Isakov and coworkers find 
the startling result that the dynamic exponent z varies 
(and apparently in a continuous fashion) along the phase 
boundary between the disordered, topological phase and 
the ordered phase. As a check they confirmed that at the 
special point where the classical eight-vertex model be- 
comes two decoupled classical Ising models, the classical 
dynamic exponent is z ~ 2.196 (which is consistent with 
the known value of z for the classical 2D Ising model). 
Thus, these classical Monte Carlo simulations imply that 
the quantum theory is unstable and that the dynamics 
of the effective field theory is not z = 2 as it is along the 
six- vertex line. 

These results motivated us to study the full quantum 
Lifshitz action and the effect of perturbations as one 



moves away from the six-vertex line. We do this by per- 
turbing directly the quantum Lifshitz model action with 
both local charge and magnetic excitations. We should 
note that this approach does not enforce the RK condi- 
tion and hence it is not equivalent to the classical sim- 
ulations of Isakov and coworkers. By an analysis of the 
RG equation we will see that the perturbed quantum Lif- 
shitz model has one line of fixed points. We find that the 
dynamic critical exponent is z — 2 along this line and 
it does not change. On the other hand, the magnetic 
excitations are described by a marginally relevant per- 
turbation. Thus this theory has two marginally relevant 
operators whose coupling constants grow monotonically. 
Here two things can happen: one option is that the quan- 
tum phase transition becomes first order, in which case 
there is no quantum criticality. The other option is that 
the system flows to a new fixed point, which is not per- 
turbatively accessible. Presumably this new fixed point 
will be have "generic" Z2 quantum phase transition with 
dynamic exponent z = 1. Clearly, our one- loop RG can- 
not be used to construct this theory. 

In the regime in which our RG calculation is reliable, 
it is consistent with a dynamical exponent z = 2. On 
the other hand, in order to reproduce the simulation 
done by Isakov and coworkers along the dual critical 
line of the classical Baxter model would require to con- 
struct a Martin-Siggia-Rose supersymmetric effective ac- 
tion which must have an exactly marginal operator to 
describe the equilibrium spatial correlations (with vary- 
ing exponents) and with a continuously changing z. Al- 
though to the best of our knowledge this theory does 
not yet exist, we can speculate on a possible scenario to 
understand the results of Isakov et al. It is well known 
from the work of KadanofP^l that the critical line of the 
Baxter model dual to the six vertex model line can be 
mapped to the partition function of a Gaussian model 
perturbed by electric and magnetic charges. The prob- 
lem that we studied in this paper is the quantum analog 
of this problem. Now, a Langevin equation constructed 
to study the classical critical dynamics of this problem 
will be local in the scalar field 4> of the 2D classical gaus- 
sian model. Since the electric charges are represented by 
local vertex operators of the Gaussian model, these exci- 
tations are also local in the stochastic process described 
by the Langevin equation. However the magnetic charges 
(or vortices), which are dual to the electric charges, are 
non-local in terms of the field <fi which represents the 
height variable. They are physically edge defects in the 
lattice picture. This suggests that the relaxation pro- 
cesses become increasingly non-local as the coupling to 
the magnetic excitations increases. Thus it is conceiv- 
able that this non-locality may give rise to a varying z. 
However, even if these considerations are correct, a con- 
tinuously varying z is inconsistent with a theory with a 
local action with an exact marginal operator. 

In the next section we provide details of the quantum 
Lifshitz model for the reader unfamiliar with the subject 
matter. In section IV, we will then return to this question 
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of flowing dynamical exponent and present details of our 
analysis. Detailed calculations can be found in Appendix 
B. 



III. QUANTUM LIFSHITZ THEORY 



Ardonne et. al generalized this construction to cases 
in which the wave function is not an equal-amplitude 
superposition as in the "RVB" wave function, Eq.(3.2|, 



but to one in which the configurations are assigned a local 
weight. In the case of the quantum eight-vertex model 
along the RK lines the wave function ia^U 



The quantum Lifshitz model is the simplest contin- 
uum field theory that has been proposed to describe the 
quantum six vertex and quantum dimer models on the 
lattice. The simplest of illustration of these ideas can be 
seen with the quantum dimer model on the square lattice 
which acts on the configuration space of a classical model 
of close packed (one dimer per vertex) dimers on a two di- 
mensional lattice. The same considerations apply to the 
quantum eight vertex model that we will discuss below, 
briefly. For a much more thorough discussion relating 
the quantum Lifshitz theory and the six vertex model, 
we refer the reader to Ref. [21]. 



\y qSv ) =^2a N ^b N "^b N ^c N ^d N ^\C) (3.3) 

{C} 

where {C} are the configurations of the eight-vertex 
model with the weights shown in Fig. [T] The norm 
(squared) of this wave function, H^gSull 2 ! is equal to the 
partition function of the classical 2D eight-vertex (Bax- 
ter) model with weights a 2 , b 2 , c 2 , and d 2 . In Figs, [l] 
and [2] we show the vertices and the resulting phase di- 
agram represented by this wave function (for the case 
a = b = 1) in terms of c 2 and d 2 . For similar expressions 
for generalized dimer models see Refe. 31 .'!2j andl40]. 



A. Quantum Dimer Models 

Let {|a)} represent a set of linearly independent and 
orthonormal states which in the dimer models are the 
dimer coverings of the lattice and in the eight vertex 
models are the allowed arrow configurations. The quan- 
tum Hamiltonians are local hermitian operators that act 
on this configuration space. The the kinetic energy term 
describes a quantum mechanical process mapping two 
configurations that differ by local moves (flips) while the 
potential energy term of the Hamiltonian penalizes each 
flippable plaquette i of the lattice. It is conventionally 
written in terms of operators acting on all the plaquettes 
of the lattice 

i 

Here Fi is the local kinetic energy term (the flip oper- 
ators) and Vi are the local potential energy terms. For 
special choices of parameters each term of this Hermi- 
tian operator obeys the RK condition^, (Vi — Fi) 2 = 
2(Vi — Fi). Hence, the Hamiltonian can be written as 
H = \ J2i Q\Qi where Qi — (V — Fi). Hamiltonians of 
this form necessarily have eigenvalues E > and more 
importantly, if one can find a state annihilated by Qi, 
then one has found the ground state. The equal ampli- 
tude sum over all possible classical dimer coverings is one 
such state. Then, the properly normalized ground state 
wave function is given by, 

iv*>> = -^Ei a >- ( 3 - 2 ) 

v a 

Z is simply number of all dimer configurations on some 
finite two dimensional lattice and is precisely the classical 
partition function of two dimensional dimers. 



B. Mapping Quantum Dimer and Vertex Models 
to the Quantum Lifshitz Model 

This analysis extends to the computation of equal 
time correlation functions in the ground state, and one 
finds that these correlators are simply given by correla- 
tion function of the classical two dimensional theoryP^ 
Hence, for T-Lrk on the square lattice, the equal time 
static correlation functions are algebraically decaying and 
%rk is said to sit between two ordered phases of the 
dimers. On the triangular lattice, however, the equal 
time correlation functions in the ground state are instead 
exponentially decaying and the quantum dimer model is 
gapped. The system is in fact in a topological phased 

The dimer configurations (as well as the six vertex con- 
figurations) can be mapped to a height variable on the 
dual lattice.^ This mapping holds for both the classical 
and the quantum problem. The dimer-height mapping is 
exact provided the dimer configurations are fully packed, 
i.e. every site of the lattice belongs to one and only one 
dimer, and that the lattice is bipartite. The violation 
of either condition is mapped in the height model into a 
topological defect which is referred to as a vortexP^ The 
same issue arises in vertex models in which the vertices 
which violate the local conservation law of the six ver- 
tex model are also represented by vortices in the (dual) 
height model picture. 

Let us assume for now that no such topological singu- 
larities are present. We can no w use a coarse-grained pic- 
ture of the (dual) height modelP^^In this effective con- 
tinuum limit, the coarse-grained height variable becomes 
a free boson tp(x,t) (see Figure [T]) with compactification 
radius given by the modulo 4 condition on the lattic e. A 
continuum Hamiltonian conjectured by Henley^^y anc j 
expanded by Moessner and collaborators^ to belong in 
the same universality class as the square lattice quantum 
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dimer model is, 

n -- 



d 2 x 



n 2 

2 



-(vV) 5 



(3.4) 



where II = p as usual. The associated Euclidean action 
is given by, 



S = 



d 3 3 



<£_ 
2 



-(vV) 5 



(3.5) 



This action also arose in the study of three dimensional 
classical statistical mechanics in the study of Lifshitz 
points in smectic liquid crystaLP^l and hence the theory 
was dubbed the quantum Lifshitz modelPf More impor- 
tantly, the ground state wave function is related to a two 
dimensional Euclidean free boson and represents a line 
of critical points parameterized by the constant 

This can be done by quantizing the Hamiltonian by 
imposing the equal time commutation relationspS 



[p{x)Mv)] = id^(x-y) 1 



(3.6) 



so that in the Schrodinger picture the canonical momen- 
tum can be written as H(x) — i6/5ip(x). Schrodinger 
equation can be written as, 



d 2 x 



6 

Sip 



$o[<p] = Et> [ 



(3;7) 

Now, the Hamiltonian can be written as a self-adjoint 
operator as in the RK Hamiltonian. Defining the opera- 
tors 

— ( — 

(3.8) 

-(-- 
V2\ Sep 



Q(x) 
QHx) 



KV 2 (p 



nV 2 p> 



The normal ordered quantum Hamiltonian is then given 
by, 



H = 



1 



d'x <QHx),Q(x)} -e vac V 



d z xQ\x)Q{x) 



(3.9) 



where V is the area and e vac is the UV divergent vacuum 
energy, 

e vac = -^]imV 2 s 5^(x-y). (3.10) 

The normal ordered Hamiltonian is a self-adjoint opera- 
tor with energy E > and therefore any state for which 
Q^oM = must be the ground state. This is simply a 
first order functional differential equation which is easily 
solved by, 



1 



e 2 



fd 2 x(V V ) 2 



(3-11) 



where Z is the normalization factor, 



Z = J [Dtp]. 



(3.12) 



Correlation functions at equal time in the ground state 



are 



(va.c\0[<p(x 1 )]...0[<p(x n )]\v&c) = 

[Dp] 0[<p(x 1 )}...0[<p(x n } e-^ d2 ^ v ^ 2 ,(3.13) 



and are hence directly related to the power law correla- 
tion functions of the free boson field theory. The behav- 
ior of the correlation functions of this theory is controlled 
by the parameter k of the quantum Lifshitz Hamiltonian. 
As in the classical 2D XY-modelpSl and in the Luttinger 
model of ID fermionic systems, there is a line of critical 
points characterized by the parameter n. The parameter 
k is related in a non-universal way to the microscopic 
models. For instance, Ardonne et. al gave an explicit 
expression for the relation between k and the weights c 
and d in the quantum eight-vertex model.^ 
In particular, the operator 



O q [p{x)} 



Jqtp(x) 



(3.14) 



for q e Z creates a bosonic coherent state which we call a 
charge excitation. This operator enters in the expression 
for the order parameters of the ordered phases of the 
quantum dimer mode^ and of the quantum eight-vertex 
modelJ^H This correlation function obeys a power-law of 
the form 



(O q (x)0. q (y)) 



1 



y 



q 2 1 4ttk, 



(3.15) 



which follows from a similar set of manipulations per- 
formed in the XY-model and using the quantities in Ap- 
pendix |Bj It follows that the scaling dimension for the 
operator O q is 

The "magnetic charge" (vortex) operator is given by^ 
6 m {x) = e i2m f d2g ar s(?-?)n( 2 -)^ ( 3 _ 17 ) 

The action of this operator on the ground state is to shift 
in the boson configuration to another with a vortex singu- 
larity at x with topological charge (vorticity) m. In other 
words, it introduces a jump discontinuity in the height 
variables and corresponds to an addition of "source" and 
"sink" vertices to the six vertex model to give the full 
eight vertex model (see Figure [l]) The defect (Dirac 
string) associated with these magnetic vortices corre- 
sponds to the non-local term exp (J dz &rg(x — z)n(z)). 
It amounts to a singular gauge transformation on the 
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field (p. Its action on the ground state wave function can 
be written as, 



{O mi (x 1 )...O mn (x n )) 



1 



Dip e 



-k f d 2 x (Vtp+A) 2 



(3.18) 



where A satisfies, 

n 

e^ViAj =2irJ2 m e S 2 (z-x e ) (3.19) 
e=i 

This expression is non-vanishing only if the charge- 
neutrality condition, m i — 0) is obeyed. This condi- 
tion reflects the fact that the wave function of the quan- 
tum Lifshitz model is vortex-free. The same considera- 
tions apply to the electric charges. 
Eq.(3.19) has the solution 



Aj = mi dj arg(z — xt) (3.20) 

t 

By using the Cauchy-Riemann condition, 

dj arg(z - Xf) = tijd 1 log \z - x e \ (3.21) 

its straight forward to find the correlation function. Spe- 
cializing to the two point function, the result is, 



(6 m (x)6_ m (y)) 



1 



V 



(3.22) 



as expected, the equal time critical exponent is given by, 
A m = 2iTKm 2 . 

From this close connection with conformal field theo- 
ries, its easy to see that the static behavior of the con- 
formal quantum critical points and specifically the quan- 
tum Lifshitz theory are well understood. The dynamic 
behavior, however, is another story and has remained 
a largely unexplored. In recent numerical computations 
based on simulations of the Langevin equation, it was 
proposed that by perturbing the Lifshitz theory by al- 
lowing for line defects (i.e. the operator O m ) it would be 
found that the dynamical exponent would flow continu- 
ously from z = 2 to z = 1. If true, this would contradict 
much of what is known. In most models the dynamical 
exponent is fou nd to jump discontinuously between inte- 
ger values PHHI However, because of the close relationship 
between the quantum Lifshitz theory and the XY-model 
where all critical exponents vary continuously, it was hy- 
pothesized that perturbations in the quantum Lifshitz 
theory would result in a flowing dynamical exponent P2l 
In the remainder of the paper, this issue is explored in 
further detail. 



IV. PERTURB ATIVE RENORMALIZATION 

In this section we will consider a system that consists of 
a 2 + 1-dimensional quantum Lifshitz model plus a set of 



local electric and magnetic perturbations. Such a generic 
perturbed system does not (in general) respect an RK 
condition and therefore it is not equivalent to a classical 
system with a partition function perturbed by similar 
operators. In other terms these two problems are not 
in general equivalent as the 2 + 1-dimcnsional theory in 
general is not consistent with the supersymmetry implied 
by the Langevin equation. 

It has been proposed that the square lattice quantum 
dimer model i s described by the pure quantum Lifshitz 
theoryP32H2U We wish to study the theory near the six 
vertex model and perturb into the full eight vertex model 
by allowing for magnetic vortices which break the ice 
rule of the six vertex model. This point is Kosterlitz- 
Thouless (KT) like and to move deeper into the ordered 
phase corresponds to varying the fugacity of the vertex 
operators. The theory along the quantum six vertex line 
is analogous to the sine-Gordon model and was studied in 
detail by GrinsteirPS where it was called an anisotropic 
sine-Gordon model. Specifically we consider one set of 
perturbations, 

(X OL 

^(O q (x,t) + 0_ g )(x,t)) = -j cos (</¥>(*,*)) (4.1) 

As noted above, the equal time scaling dimension of this 
operator is A q = q 2 /(87m). 

Secondly, we need to perturb the quantum six vertex 
model by allowing the ice rule to be broken. This corre- 
sponds to an operator, 

^(d m (x,t) + G- m (x,tj) (4.2) 

where the time dependent vortex operator is given by, 
6 m {x l ,t l )= (4.3) 
exp (^i2m J dr J d 2 z arg(z — xc)8{t — ti)Tl{z, r) 

This operator gives the same equal time correlation func- 
tions (3.22) and its equal time scaling dimension is given 



by A,„ = 2-nnm 2 . Analogous to the anisotropic sine- 
Gordon modelj^ the coupling constant at the KT-point 
is k c = l/(327r), and one easily finds that the O q is 
marginal for q = 1 and O m is marginal for m — 8. 

Taking the effective low energy theory to be the quan- 
tum Lifshitz theory, we consider the action, 



S = 



d 2 x 



<lt \{d t p) 2 



-iy 2 ^ 



A(Vv?) 
a 
2a 



4 (<V 



■ <5(VV) 2 + u(V<p) 1 
a 

2a~ 4 



,4" 



(4.4) 



We expand about the fixed point where both O m and O q 
are marginal at k c = 1/327T. Defining 8 to be the distance 
away from marginality so that 6 <C 1, we look at the 
quadruple expansion in the small parameters S, u, a and 
a while treating A = A(u, a, a, 5) as a renormalization 



group parameter; one fixes A = to give the Lifshitz 
theory. In this sense it is a parameter that constrains the 
flow equations properly as opposed to a coupling constant 
that flows under the RG. While u is marginally irrelevant, 
the justification for including such an operator is that it 
is generated under renormalization anyway. 



A. Renormalization group for the perturbed 
quantum Lifshitz model 

Now, we want to study the flow of the various coupling 
constants when the variables u, a, a and 5 are small. This 
is done through a perturbative renormalization scheme. 
In such a procedure, one begins with an effective low en- 
ergy field theory described by an action So and a set of 
operators Oi with coupling constant gi and scaling di- 
mension Ai that will be treated perturbativcly. Treating 
the coupling constants as small, one can expand in pow- 
ers of <7i, 

Z = Tr e -^-S^ s T,i9iOi{S) 
~Z 



1-1 d 3 a? 

i 

^ J d 3 xd 3 yJ2g l gAo t (x)OAy)) 



(4.5) 



where the correlation functions are taken with respect to 
the low energy effective field theory. The renormalization 
is then performed by rescaling the short distance cutoff 
a — > Xa where A = (1 + 5) where S > and asking 
how the couplings rescale so that the partition function 
is preserved. In the first term, this rescaling is simple. 
The coupling constant g t rescales as, 



Hi 



x d - 



'9i 



g t + (d- Ai)gi8 



(4.6) 



In the second, one may use the operator product expan- 
sion. Under a rescaling, the limits of integration can be 
written as 



/ -/ 

J\x-y\>(l+8)a J\x-v\: 



(4.7) 



\x—y\>a J a<\x— y\<a(l+S) 



The first term only gives the original contribution to the 
Hamiltonian. In the second, one may apply the operator 
product expansion. We know that when two operators 
are nearby, we know that they can be thought of as ef- 
fectively "fusing" into another operator in the theory. 



x + y 



(4.8) 



where C\ v are the fusion coefficients. 

To find the fusion coefficients, one looks at the three 
point function. Its a generic feature of d-dimensional crit- 
ical systems that the three point correlation functions are 



restricted by translation, rotational and scale invariance 
to have the form 

(Oiix^O^Okixa)) ~ 



1 



(Ai 2 ) Vi+y j~ Vk (Ai 3 )^ i+ ^ fc ~» (A 2 3) yj+yk ~ Vi ' 



(4.9) 



where = \xi — Xj\ and yi is the scaling dimension of 
Oi. The added complication in anisotropic fixed points 
is that scaling invariant functions of s — \x\ /t may also 
appear, but the structure remains the same. Defining 
Sij = Afj/tij, generically what we find that the three 
point function are scaling functions of the form, 

(0^)0^X2)0^X3)) ~ 



/(S12)/(S1 3 )/(S 23 ) 



(Al2) 



Vi+Vj-Vk 



'(A 



lz )Vi+Vh-Vj (A 2 3) w+2 "=~ 



(4.10) 



Now, in the limit that Ai 2 <C A13 and A i2 <C A 23 and 
analogously for tij, the three point correlation function 
factors. If 27? = X1+X2 and 2f— x\~xi (and likewise for 
time) while further setting the dummy variable X3 = 0, 
this can be written as, 



(O,(f 1 )O i (f 2 )O fc (0)) ~ c ijk 



f(s r ) P(S R ) 
r Vi+Vj-Vk J$Vk 



(4.11) 



where the cyfc are numbers and the functions f(s r ) are 
scale invariant functions of s r = r 2 /t. Because of the 
more complicated scaling functions present, less diver- 
gent terms in R are also present, but we keep only the 
leading divergent behavior here. To connect with the re- 
sults in GrinsteirPS, we finally want to integrate all values 
of s to some fixed value sq. Then, by taking the time to 
be small at this point, one integrates < s < 00. In ad- 
dition, exchanging the time variable for a dimensionless 
s variable, one needs to only consider terms which are 
divergent as r~ 4 in the limit r — > 0. 

By doing so, one finds that the contribution coming 
from the terms second order in the coupling constants, 
yields a contribution, 



IE 



,Ai+A,/-A,- 



ds r dClf(s r ) 



d 6 xd 3 y ( O 



a<\x-y\<a(l+S) 



.(4.12) 



As mentioned previously one integrates over all possible 
values of the aspect ratio and the solid angles. Provided 
that the integrals converge, this is a legitimate proce- 
dure and yields a simple number. Re-exponentiating the 
result, one finds that the coupling constant then flows as, 



dgj 
dS 



1 



(d - Aj) gj - ^S d 2^ e itil gigi, 



<■[, 



(4.13) 
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Hence, to find the /3-functions various operator fusion 
coefficients cr? must be computed. This is outlined in 



Appendix B. Defining the quantity, 

s/4 



Ji 



ds s 



-2 + 2i 



exp 



-dx 



(4.14) 



so that numerically, J\ = 0.290 and J 2 = 0.111, one 
finds that for the quantum Lifshitz model, the non-zero 
operator product expansion coefficients are, 



(V^) 4 (f) :: (V^) 4 (y) 





4 


' 2 








log 




) a 4 




12 5 





5 , /4 
-12 +l0g 3 



51 : (VV) 2 (^) 



« ,^ 1 e 4 ^ J 2 <W 

4 4' 7r a 4 

P 4 7 



e 4 ~< J 2 6£ .„ 



1 e 4 ^ J 2 51 



8 (f)0_ 8 (^~ 7 _^_:l 



4 4 7 



e 4 ^J 2 « 2 - 

4^32^ :(V ^ } {R) 



4 7 



Ji« : (V 2 ip) 2 (R) 



O x (x) : (V^) 4 (y) 
1 (x) : (VV) 2 (y) 
6 8 (£) : (V<p)\jf) 

o 8 (f) : (vV) 2 (y) 



46 

1 16, /32 
2^V l0g U 



1_. 

32tt C 

1 16 
2^V l0g V27 



: O-i(iZ) 
32 s 



: O-i(ii) 



32tt 



: 0- 8 (i?) 



a z 1D7T 



(V^) 4 (i?) : 



e 4 7 
46 



ttJ 2 <« : (V<p) 4 (i?) : 



e 4 7 
46 



7rJ 2 <5£ : (V<p) 4 (i?) : 



(4.15) 



From the operator product expansion and fusion coef- 
ficients, one can easily find the /3-functions for the various 
coupling constants. The renormalization group equations 
are similar to t hose of the two dimensional isotropic sine- 
Gordon modelPSHSHZl The /3-functions have the form, 

4 7 



£(«) 



41 



32tt 

/3(a) = -4Sa - 



16 



7re 
log 



J2 (a 2 



a 2 ) 



(3(a) = 4(5a + — log 



-Jl 



(a 2 - a 2 ) 



log 



5 

12 



1 

47 •"■ V" -1 2^2 ^ 

(4.16) 

Ignoring (3(a) the anisotropic sine-Gordon model was 
studied to one loop order by GrinsteirpSI and our re- 
sults qualitatively agree with his. There are small dif- 
ferences in the numerical prefactors, but we suspect that 
this arises from differences in the renormalization proce- 
dure. It is not clear that the approximations used here 



amount to the same approximations used by GrinsteinPHl 
As mentioned previously, A — A(u,a,a) is a renormal- 
ization group parameter and determined by the condition 
that ,9(A) = 0. 



A(u, a.a) r - 



log 



2 

12 



4 5 32 



a 2 ) 



(4.17) 



In particular the fix point of the theory is still the trivial 
fix point, (S*,u*,a*,a*) = (0,0,0 0). 

The one-loop RG equations, eq.( |4.15 1, have the follow- 
ing structure. They have one critical line at a = a = 
and u = 0. This is the six-vertex model line. Along this 
line the system is represented by the quantum Lifshitz 
model which has an exactly marginal operator, (y 2 ip) 2 
(and d t ip) 2 )- Along this entire line the coupling constant 
u — > very slowly (since (Vy;) 4 is a marginally irrelevant 
operator.) On the other hand, near the end of the six- 
vertex line the magnetic excitations with coupling con- 
stant a > 0, are also marginally relevant, and this cou- 



10 



pling constant grows under the RG. Thus, instead of a 
second hxed line, the phase boundary for a > is either 
a first order transition (in which case there is no quantum 
criticality) or it flows to another quantum critical point 
which is not perturbatively accessible. Our one-loop RG 
cannot solve this problem. Since the z = 2 dynamics is 
special and "non-generic" one possible solution is that 
the quantum critical theory is controlled by a z ~ 1 fixed 
point (possibly with Ising criticality). 

Treating A as a renormalization group constant and 
fixing it to be A = for the z — 2 quantum Lifshitz 
theory, it is easily seen that the renormalization group 
equations share many of the same properties as the XY- 
model. A critical surface is described by the equation, 
fu 2 = g(a 2 — a 2 ). In the limit where there are no mag- 
netic vortex excitations, the /^-functions describe a line of 
fixed points and the charged excitations are marginally 
irrelevant operators. This is in agreement with earlier 
studies of the anisotropic sine- Gordon modelP With 
the charged excitations suppressed, another line of fixed 
points is described by the /3-functions and the magnetic 
vortices are marginally relevant operator that drive the 
system away from the critical poin t, much like in the 
generalized sine-Gordon modelpSB3 Like the generalized 
sine-Gordon model, this has two regimes depending on 
the sign of S. For 6 — > — oo, the vertex operator is rel- 
evant and the vortex operator is irrelevant: a grows to 
large values while a tends to zero. In this case, the system 
is locked into the phase where for a < 0, <p = (2n + l)ir 
while for a > 0, ip = 0. In the opposite limit, 5 — > +oo, 
the opposite is true. Here the vertex operator is irrelevant 
and vortex operator is relevant. Now it is a that tends to 
zero while a grows large. Unfortunately, the argument is 
more complicated and the field IT can not be thought of 
as being simply locked into a specific configuration. 



B. RG flows and quantum stability analysis 

Finally, one can linearize the /3-functions about the 
critical point. The renormalization group equations have 
the structure, 



% 
dr 



Cijkgjgk, 

j<k=i 



(4.18) 



where <?i = u,g 2 — a,g 3 = a and g 4 = S and t = 
log|x|/a parameterizes the scale factor by which one 
changes the short distance cutoff. Now, suppose that 
one is close to the critical point and that a, a <C u. 
Then at least for small r near the critical surface, one 
can solve the renormalization group equations by elim- 
inating terms 0(a 2 ) and 0(a 2 ). This leaves the set of 
two coupled differential equations, 

du n 2 



Ot 



which have the solution u = w (l + ChiUqt)^ 1 and 
5 = 5q + C^iiC-f^u where uq,Sq are arbitrary integra- 
tion constants. It is short work to solve the remaining 
differential equations, 



da 

da 
cV 



(-4(5 - C 2 i2u)a 
{AS + C 313 u)a, 



from which one readily finds, 



a = ao(l + C nl u Q T) x e- 4S( > T 
a = a {l + C nl u T) y e 4S <> T , 



(4.20) 



(4.21) 



-C411U , 



(4.19) 



where x = — {ACm + C212C111) /C m ~ —1.96 and y = 
—x. 

For the moment, we restrict ourselves to the region of 
parameter space where a — 0. With 6 = 0, the fixed 
point u = 5 = a = fixed point is reached as the scale 
r — > 00. Specifically, as in agreement with GrinsteinpHl 
we find that asymptotically, u ~ A u /t,S ~ Ag/r and 
a ~ A a r x where A u ~ C^ n ,As ~ C4nC m and Ag ~ 
oto(CmUo) x ■ For S < one finds that a flows to infinity 
as the scale is changed. Flows stay close to the critical 
surface when r is not too large. To be precise, for r ~ 
I /So the flows will remain near the critical surface as 
a does not flow to infinity too quickly. Recalling that 
r = log |af|/a this gives an idea of how large the critical 
domains are in this phase and an idea of the correlation 
length. One finds that, £/a ~ exp [f(5)~ a ] where the 
critical exponent a — 1. On the surface where a = 0, 
similar considerations leads to identical conclusions only 
the regions S > and S < are exchanged. 

Finally, we examine the situation where a, a 7^ 0. 
There is no critical point here since as r — > 00 either 
a flows to larger values (when 6 < 0) or its dual a flows 
to larger values (when 6 > 0). However, again for val- 
ues of r that are modest, the flows stay near the critical 
surface. Focusing on the region S > 0, one finds that for 
r ~ j^, a does not flow too quickly away. Once again, 
this gives a rough estimate to how the correlation length 
scales, namely, we find again that £/a ~ exp [/(i5) _tT ] 
with cr = l again. Similar arguments may be applied to 
6 < 0. For more details, see Figure [3] 

The main conclusion of the RG analysis is that the per- 
turbed quantum multi-critical point of the quantum Lif- 
shitz model is that when magnetic excitations are present 
they lead to a runaway behavior of the RG flows. Since 
at this point the theory has two marginally relevant per- 
turbations the runaway process is very slow inducing log- 
arithmic corrections-to-scaling in the correlators. The re- 
sult of the runaway flow is that the quantum phase tran- 
sition between the ordered phase and the Z 2 topological 
phase is either a (fluctuation-induced) first order transi- 
tion, or a continuous quantum phase transition controlled 
by a fixed point (presumably with 'generic' z = 1 dynam- 
ics) which is not accessible to the perturbative renormal- 
ization group theory that we have presented here. This 
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FIG. 3. The renormalization group flows in the full quan- 
tum eight vertex model. The renormalization group flows are 
reminiscent of Kosterlitz-Thouless physics. Dashed lines are 
flow lines in the plane where a = while the solid lines are 
in the plane with a — 0. There is a single fixed point at 
(S,a,a) = (0,0,0). 



result is in markedly different from the conclusions of Ar- 
donne and coworkers,^ based on the equal-time behav- 
ior of the correlators in the Baxter wave function, and 
of the classical Monte Carlo simulations of fsakov and 
coworkers In the next section we analyze these issues. 



V. DYNAMICS AND THE RK CONDITION 

In the previous section, we examined the parameter 
space of the z — 2 quantum Lifshitz model in the pres- 
ence of both magnetic and electric charge perturbations 
to one loop order. There we found that the magnetic 
perturbations at the quantum Baxter multicritical point 
lead to a runaway RG flow of the coupling constants 
which render the quantum phase transition either weakly 
(fluctuation-induced) first order or continuous but con- 
trolled by a non-perturbative fixed point. This result 
is markedly diffe rent tha n what happens in the classi- 
cal Baxter model F' 46 ' 47 ! We will see that the quantum 
dynamics is different deepening on whether the RK con- 
dition is preserved or not 

Recently, there has been renewed interest in this 
model. Using classical simulations Isakov and cowork- 
ers found that the dynamics of the quantum Lifshitz 
model along the phase boundary between the ordered 
and topological phases of the model was in fact unsta- 
ble and flowed from z = 2 to values z > 2 (albeit in a 
non-monotonic way) The simulations along this phase 



boundary are classical in the sense that they follow the re- 
laxational dynamics of the classical Baxter model. How- 
ever, in the sense that we will now describe, these simu- 
lations nevertheless represent the quantum dynamics of 
the quantum eight vertex model of Ardonne et aZP^ 

The simulations of Rcf. [29] use the fact that Hamiltoni- 
ans that obey the RK condition can also be used to repre- 
sent classical stochastic processes of relaxation to thermal 
equilibrium. For any discrete classical model, the classi- 
cal dynamics can be described by a master equation.SsHini 
If p a is the probability of having a specific configuration 
of dimer (or arrows in the eight vertex model) a, then 
the dynamics can be written down as a Master Equation 



d t p a {t) = ^W aP pf}(t) 



(5.1) 



where the sum runs over configurations (3 that differ from 
the configuration a by a single dimer flip. For the classi- 
cal dim er mo del, the (positive) transition matrix W a p is 
given bjEUm 



a/3 



V a 8 a p - F ( 



a/3, 



(5.2) 



where V a counts the number of flippable dimers there are 
in a configuration a and F a p dimer flip operator. 

The connection between the classical dynamics repre- 



sented by the Master Equation, Eq.(5.1) and a quantum 



Hamiltonian (acting on the same configuration space) is 
as follows.^ Let p a (t) be the classical probability of a 
configuration a at time t, and let w(a) be a classical 
equilibrium Gibbs weights for the configurations a. It 
then follows that, provided the transition matrix W a ,p 
obeys the detailed-balance condition, 



(5.3) 



(no sum over repeated labels is implied) , then in the long 
time limit the system is guaranteed to reach the thermo- 
dynamic equilibrium state , 



lim p a (t) = w a 

i-»+oo 



(5.4) 



Now, upon a simple rescaling of the transition matrix, 

W a , P = W^Wa^W 1 / 2 (5.5) 

and of the time-dependent probabilities, 



Pa(t) 



w. 



-1/2 



Pa(t) 



(5.6) 



allows us to identify the transformed Master Equation 
with a quantum mechanical evolution in imaginary time 
with a (hermitian) Hamiltonian 



H a ,p = ~W a ,p 
and the wave functions become 



(5.7) 
(5.8) 
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In particular the ground state wave function of this 
Hamiltonian has the form of a sum of local Gibbs weights, 

l*o>=E< /2 |«) ( 5 - 9 ) 

a 

By construction this state has exactly zero energy, 
H\^ ) = 0, and the norm (squared) of the wave func- 
tion is the classical equilibrium partition function, 

ll*o|| 2 =5>« (5.10) 

a 

It is easy to see that the constraints required for the clas- 
sical stochastic process to reach the equilibrium Gibbs 
state are equivalent to the RK condition that the Hamil- 
tonian is a sum of projection operators. This mapping 
also implies that the spectrum of relaxation times of the 
classical evolution problem is the same as the spectrum 
of energy gaps in the quantum problems (with the gaps 
being the reciprocal of the relaxation times). This re- 
lation between classical stochastic processes and quan- 
tum Hamiltonians that satisfy the RK condition has been 

exten sively used to construct non-trivial 2D quantum 
sta tes.™mil3] 

Isakov and coworkers^! used the relation between clas- 
sical relaxation processes described by a Master Equation 
and 2D RK quantum Hamiltonians to carry out extensive 
simulations. Here we will focus on their results for the 
quantum eight vertex (Baxter) model of Ref. 21 which 
obeys the RK condition. 

Along the six vertex model line, where the arrow con- 
figurations that violate the local conservation law are not 
allowed (hence the weight d = 0) the classical simulations 
find that the dynamical exponent is z — 2 as predicted 
by the quantum Lifshitz model.^ This result also agrees 
with our findings in the preceding section in the absence 
of magnetic excitations, w hich rep roduce the earlier RG 
results in the same regimeP32SE2l 

On the other hand, along the phase boundary be- 
tween the ordered phase and the Z2 topological phase, 
parametrized by c 2 = d 2 + 2, their classical simulation 
obtains a dynamical exponent z > 2. This result is con- 
sistent with the fact that the classical dynamics along the 
six vertex model line, which has a local conservation law 
and consequently the quantum Hamiltonian (and hence 
the Liouville operator) has a continuous symmetry. In 
this case the classical critical dynamics is in the model 
B class in the Hohenberg-Halperin classification,^ and 
hence has z = 2 along this entire line. In contrast, along 
the phase boundary between the Z2 topological phase 
and the the ordered phase the local conservation law is 
broken but the quantum Hamiltonian has a local Z2 sym- 
metry. Thus in this case the universality class of the clas- 
sical critical dynamics is expected to be model A. Their 
result of z > 2 is consistent with the inequalities that 
the dynamical exp onent is known to obey in the case of 
Model A dynamics! 45 * 52 ! Moreover at a special point on 
this phase boundary, (c 2 ,d 2 ) = (y/2 + 1,^/2— 1), the 



norm squared of the ground state wave function is equal 
to the partition function of two decoupled 2D classical 
critical Ising models. The value of the dynamical critical 
exponent obtained by Isakov and coworkers turns out to 
be z — 2.196 which is consistent with the best estimates 
for z in the 2D critical Ising model with model A dy- 
namics. However what is surprising is the fact that they 
find an exponent z which varies continuously (and non- 
monotonically) along this entire phase boundary, rising 
from a value close to 2 at (c 2 , d 2 ) = (2, 0) (the KT tran- 
sition of the six vertex model) , peaking at the decoupled 
Ising point mentioned above and slowly falling off past 
that point to values still obeying z > 2. 

Isakov and coworkers also considered the effects of per- 
turbations to the eight vertex model Hamiltonian (at the 
Ising decoupling point) that violate the RK condition. In 
this case they did quantum Monte carlo simulations and 
found that for all values of this RK-violating perturbation 
the dynamics obeys z = 1, instead of z — 2.196 when the 
RK condition holds. The value z = 1 is consistent with 
the critical behavior of the 2D quantum Ising model (in 
a transverse field), whose quantum Hamiltonian is (es- 
sentially) the same a transfer matrix of the 3D classical 
Ising model. The quantum Ising model is well known to 
be the dual of the 2 + 1-dimcnsional Z2 (Ising) gauge the- 
ory. This is also consistent with the fact that the Kitaev 
Toric Code point is equivalent to the ultra-deconfined 
limit of the Z2 gauge theory. This result suggests that 
perturbations to the 2D quantum Hamiltonian that vi- 
olate the RK condition lead to a system with 'generic' 
z = 1 'relativistic' dynamics. 

This result of Isakov and coworkers is consistent with 
the renormalization group analysis presented in the pre- 
ceding section. There we showed that by perturbing the 
quantum Lifshitz model in the vicinity of the KT tran- 
sition, where both the electron and magnetic perturba- 
tions are marginally relevant and lead to a runaway RG 
flow. Although the quantum Lifshitz model itself is con- 
sistent with the RK condition, generic perturbations are 
not consistent with this constraint ,and in general kill the 
quantum criticality, although it is still possible to reach 
this multicritical point by fine-tuning of the parameters. 
However, in the presence of magnetic perturbations this 
is no longer possible as their effective coupling constants 
flow to strong coupling. Thus, in this case the physics is 
controlled by a fixed point which not accessible in per- 
turbation theory. 

Although the quantum Monte Carlo simulations of 
Isakov and coworkers were done only near the special 
decoupled Ising point, combined with our RG analysis 
they suggest that in general, if the quantum phase tran- 
sition between the Z2 topological phase and the ordered 
phase remains continuous, then it is always in the uni- 
versality class of the deconfinement quantum phase tran- 
sition of the 2 + 1-dimensional Ising gauge theory (which 
has z — 1). Contrary to the case of the quantum Lifshitz 
model, the quantum phase transition of the Z 2 gauge the- 
ory (or, equivalently the 3D classical Ising model) does 
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not have any marginal operators and hence the critical 
exponents are fixed and no longer can vary as functions 
of any parameters. 

What remains to be understood is the result of Isakov 
and coworkers that at critical dynamics of the quantum 
Baxter model which satisfies the RK condition has a dy- 
namical exponent z that can vary smoothly along the 
phase boundary between the topological and the ordered 
phases. This result is puzzling for the following rea- 
sons. Let us consider first what happens along the six 
vertex model line where the simulations obtain a dynam- 
ical exponent z = 2, which is consistent with the quan- 
tum Lifshitz scaling. The classical critical dynamics of 
the coarse-grained height field <p(x,t) can be expressed 
in terms of a Langevin equation of the forrrP^ 



(5.11) 



where the classical equilibrium state has a Gibbs distri- 
bution with weight exp(— S 1 [<£>]), and is a Gaussian dis- 
tributed noise with the properties, 

(£(*,*)> =0, 
(£(f, £)£(£', 0) =tt5 2 (x-x')8(t-i?), (5.12) 

where is the standard deviation of the noise. Along the 
six- vertex line we can take £>6t>[^] to have the (Euclidean) 
sine-Gordon form 



g n cos nip 



(5.13) 



Here g n is the coupling constant for electric charge exci- 
tations, labelled by n, in the classical problem. 

It is possible to express this stochastic process in terms 
of a path integral in d — 2 space dimensions plus time 
("Monte Carlo time" )P3 The condition that the process 
converges to the equilibrium distribution and satisfies the 
Langevin equation lea ds to a path-integral with a local 
supersymmetric action! 36 ! 38 ! In particular, it is the su- 
persymmetry that guarantees that the process relaxes to 
the chosen equilibrium state. Thus, the RK condition 
for the quantum Hamiltonian, that requires the ground 
state to have exactly zero energy, implies a form of su- 
persymmetry. On the other hand perturbations to the 
quantum Hamiltonian that break the RK condition also 
break the supersymmetry and, hence, are not consistent 
with a Langevin relaxation processes. Hence generic per- 
turbations to the RK Hamiltonian cannot be simulated 
by classical relaxation processes. 

One important property of the Langevin process of 
Eq.(5.11) for the equilibrium state defined by Eq. (5.131, 



is that it is local. In other words, electric charge per- 
turbations are local. We know from the simulations of 
Isakov and coworkers that in this case the dynamics has 
z = 2, which is the same as in the quantum Lifshitz 
model. One the other hand if we want the equilibrium 
state of the classical baxter model we need to involve 
both electric charge and magnetic charge excitations! 3 ^ 



The classical (Euclidean) action now has the form (for 
suitable choices of n and to) 



S 8v [if] = J d 2 



K 



(V</?) 2 + g n cos nip + g m cos mip 



_ (5.14) 

where to are the magnetic charges. The field <p is the 
dual field and it is related to the coarse-grained field ip 
by the Cauchy-Riemann relation 



(5.15) 



where is the second rank Levi-Civita tensor. 

The solution of the Cauchy-Riemann relation implies 
that the fields ip and its dual <p are non-local with respect 
to each other. This means that the Langevin process 
treats both types of excitations in a very different way 
and, in particular, the presence of magnetic excitations 
makes the process non-local. Thus, while the equilibrium 
partition function is invariant under duality (which ex- 
changes electric and magnetic charges) the dynamics is 
not. It is reasonable to speculate that this difference and 
the inherent non-locality in the dynamics could possibly 
lead to a value of z > 2. Nevertheless, this still does 
not explain in any obvious way why isakov and cowork- 
ers have a value of z that varies continuously along the 
phase boundary of the eight vertex model. 



VI. CONCLUSIONS 

Motivated by the work by Isakov et. alP^ and their 
surprising observation of a continuing varying dynami- 
cal exponent in the quantum Lifshitz theory, we have 
performed a perturbative RG analysis of the 2 + 1- 
dimensional quantum field theory. We added both elec- 
tric perturbations and magnetic perturbations what are 
both marginal at the same point along the six vertex line 
of the quantum Lifshitz model. In order to do the pertur- 
bative RG analysis, we have extended the operator prod- 
uct expansion to the case of anisotropic critical points. In 
the case where magnetic vortices are ignored, our results 
are in agreement with those found by Grinstein. 30 How- 
ever, some minor differences that were found presumably 
arise from differences in how the renormalization scheme 
was implemented. 

By an analysis of the RG equations we saw that the 
perturbed quantum Lifshitz theory has a single fixed 
line and that the dynamical critical exponents is z = 
2. On the other hand, the magnetic vortices were 
marginally relevant and hence the theory has potentially 
two marginally relevant operators. In such situations, 
two things may happen: (1) there is no quantum critical- 
ity and the quantum phase transition become first order, 
or (2) the system flows to a new fixed point, which is not 
perturbatively accessible. We guess that this may be the 
z = 1 theory, but our perturbative RG calculation at one 
lop can not be used to access this theory. Even given 
these caveats, our results do not seem to be compatible 
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with a continuously varying dynamical exponent as im- 
plied by the numerical simulations carried out by Isakov 
and coworkersPSl It remains to be understood how the 
how the continuously varying dynamical exponent found 
by Isakov and coworkers can be consistent with a field 
theoretical perspective. 
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Appendix A: Useful Quantities 

The Green function is given by 

G(f,t) = (<p(f,t)<p(0)) 
Defining s = \f\ 2 /n\t\ and f(s) = (1 - e" 8 / 4 ), 
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8ttk 



log 
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V 4 G(r,t) 
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Appendix B: Computation of the fusion coefficients 

In this appendix, we give a detailed account of the computation of the non-zero fusion coefficients in the OPE. 
Defining x — y = 2r and x + y = 2R and t x — t y — 2t r and t x + t y — 2tji recall that we will be interested in the 
limit where R ^> r and tn ^> t r . We are interested in the short distant (r — > 0,t r —> with s = r 2 /n\t\ fixed) 
divergent behavior of the three point functions. Useful quantities used in this section are defined as, R± = (R— z)±r, 
T± = {t R -t z )±t r and s ± = |i?±| 2 / K |T ± |. 



In this section, we provide details of the OPE for two vertex operators (0\0-i). These correlation functions can be 
computed in much the same spirit as the two dimensional, critical XY model at the KT transition. Here, we look at 
the fusion of two vertex operators into (a) the (V<p) 2 operator, (b) the iy 2 ip) 2 operator and (c) the (V<p) 4 operator. 



a. A. (Vtpf 

We point split the operator (Vtp) 2 (z) — \im Zi ^ z V Zl ip(zi)V z . 2 Lp(z2) ■ It can be shown that upon subtracting self 
energy contributions that the only non-vanishing contribution is given by, 



1 



-q G rcg (x-y) 



lim.dHa*' (22) f[ [G(x - Zi ) ~ G(y - 
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With the aforementioned variables R±,T± and s±, it is easily seen using the results from Appendix |A"| that, 

(47tk) 2 \ R\ R 2 _ R\R 2 



(di(G(x - Zi) ~ G(y - Zl ))f = j—^ ( J -^> + J -^> - +^p/( s+ )/( s _) 



In the limit where R 3> r, t R 3> t r we have the simplification s + ~ s_ — sji_ z and R + ~ i?_ = R — z. Thus arriving 
at the result, 



<0i(ar,t«)0-i(i7>*v) : (V^(z) :) ~ — ^ (fl - s) 4 (B1) 

b. Fusion into (V 2 ip) 2 

Once again point splitting the operator (V 2 </j) 2 (z) = hm 2 ._j. 2 V 2 i <^(zi)V 2 2 </?(z 2 ) and subtracting divergent self- 
energy terms, one arrives at, 

(O^x, t x )0^(y, t y ) : (VV) 2 M ^ -^e"^- «V 2 2l V 2 2 f[ [G(x - *) - G(y - *)] . 

i=l 

The non-trivial quantity to examine is then (in the limit R 3> r and t R 3> i r ), 

V*[G(ar-z)-G(|, - 



(8™) VI2 + I 

The difference in IT+I" 1 and | r_ | ^ 1 can be studied by writing \T+\ = y/(t R - t z + t r ) 2 and \T_\ = y/(t R -t z - t r ) 2 . 
Hence, we arrive at the result, 

-q 2 G rcg (r) t 2 -s R -z/2 

{O x {x,t x )0- X {y,t v ) : (VV) 2 (*) :) * r -j. (B2) 

(8ttk) 2 - i z ) 4 

c. Fusion into (V(/?) 4 

The computation of this fusion coefficient is done as the previous two cases. We point split the operator and 
subtracting the self-energy terms, its possible to show that the non-vanishing contribution is given by, 



By a similar set of manipulations, in the limit R^> r and (j;> t r , its possible to show that 



(dfoOO-iGT,*,) : (V^) 4 (z) :) ~ g — ^ ( fl _ ■ ( B3 ) 

2. Fusion of : (Vcp) 4 :: (V</?) 4 : 

In this section, details are provided for the fusion of two operators (: (V^) 4 :: (Vt^) 4 :). The non-zero fusion 
channels are (a) the (V^) 2 operators, (b) the (V 2 (p) 2 operator and finally (c) the (V<p) 4 operator. Useful quantities 

defined in this appendix are D\ = (d\ Xl) d k '^ G{x x - yi)) and dP = (d^^V^Gixi - z 2 ))- 

a. Fusion into (Vy) 2 

Point splitting each of the operators present, this channel can be computed as, 

Km, Km, Km (d^d^d^&'^A (d^d^d^d 1 ^) (d^ET^A 
x(tp(x 1 )(p(x2)tp(x 3 )(p(x 4 )^(y 1 )tp(y 2 )<p(y3)tp(y4 : )<p(zi)<f{z2)) 
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Simple combinatorics reveals that there are two classes of possible contractions for the matrices where indices are 
contracted as Tr(DD)Tr(DDD) and one where the indices are all contracted together Tr(DDDDD). The result 
follows. 



(: (V^) 4 (f) :: (V^) 4 (y) :: (V^) 2 (z) :) 



(7r«;) 2 r 3 



(: (V(^) 2 (i?) :: (V(^) 2 (z) :) 



(B4) 



72 

(*2) 



b. Fusion into (V 2 ^) 2 

As in the previous two cases, one needs to compute the quantity, 

lim lim lim f d\ x ^d^df ^d j '^A ( d { k yi) d k ^d { e V3) d e ^) ( V 2 , VI 

Xi^rx yi^ry Zi^rz V J / V / \ ^ 1> ^ 

x(ip(x 1 )ip(x2)^(x 3 )ip(x i )ip{yi)ip(y2)v(y3)^(y4)v(zi)ip{z 2 )) 
Again, one has to treat the different contraction of indices. One representative from this set is, 

(df Cl) d k ^G(x 1 -y 1 )) (d i ^d^ ) G(x 2 -y 2 )) (df 3 V<^G(x 3 - y 3 )) 
x (d^V\ Zi) G{ Xi z 2 )) (d^Vf Z2) G(y 3 z 3 )) 

There are of course many diagrams that yield the same contraction of indices. These are similar to the previous case. 
One finds that the possible contractions are given by Tr (D(r)-D(r)) xd(R + )-D(r)-d(R-) and by d(R + )-T>-G-'D-d(R). 
Hence one finds that there is a non-vanishing contribution from this correlation function to this fusion channel. 



(: (V<p) 4 (x) :: (V<^) 4 (y) :: (V^) 2 (z) :) 



{-KK) 2 r 



12 



log 



(: (V» 2 (i?) :: (V 2 <p) 2 (z) :). 



(B5) 



c. Fusion into (V^) 4 

To compute this correlation function, one needs to examine, 

Jim_ lm lim_ (d^d^d^d 3 ^) (d { k yi) d k ' M d { e V3) d e ' {Vi) ) (d^d^^d^d^A 

X(^(xi)(^(x 2 )^(x3)^(x4)^(yi)^(2/2)^(j/3)^(y4)<y5(zi)^(22)^(2:3)^(24)) 

One has to treat the different contraction of indices carefully. The combinatorics for involved is more complicated than 
the other two, but careful counting reveals that there are four channels. These are simply given by the different ways 
of grouping a contraction of six matrices: [Tr(DD)] 3 , Tr(DDDD)Tr(DD), Tr(DDDDDD) and Tr(DDD)Tr(DDD). 
In the limit R^> r, we find a non-vanishing fusion coefficient, 



(: (V#(f) :: (V^) 4 (y) :: (V^) 4 (z) :) 



41 1 



167TK r 



(: (V^) 4 (i?) :: (V^) 4 (z) :). 



(B6) 



3. Fusion of O s O- 8 

In this section, we compute the fusion of two dislocation operators OgOs- Non-zero fusion channels are in the (a) 
(V<^) 2 channel, (b) the (V 2 </?) 2 and the (c) (Vtp) 4 channel. 



a. Fusion into (Vip) 2 

Point splitting the operator (Vip) 2 , this quantity can be computed by taking functional derivatives of the quantity, 



lim d it{zi) di 2) (d 8 (x)d- 8 (y) : V (li)^ 2 )> = 



DipDTl expl / d 3 x ^-ILip 



1. 



-IT 



2,_\2 



+i4(f , i)n(f, t) + J(f, t')ip{x!, t') 
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where 

A{x', t') = 2im (axg(x' - x)S(t' - t x ) - arg(a r ' - y)S(t' - t y jj 
Integrating out II and then, redefining ip = ip' + a(x' ,t') where 

a(f , t') = 2im (axg(a? - x)Q{t' - t x ) - arg(a? - y)Q{t' - t y )j , 

one is able to eliminate the field A(x', t') at the expense of adding the field a(x' ,t'). 
Kmd it{zi) dl Z2) (6 8 (x)6- S (y) : tpfrMZ)) 

-2-nm 2 K 2 G rag (x-y) ( r i 



(B7) 
(B8) 



— cxp ^ y dV<*V - J(x')G(x' - 2/0^(2/') + ^ a {x')G(x' - y')J(y) - J(x')a(x') j 

Taking functional derivatives, there are many terms, however the only terms that survive are those were the points 
Zi are not connected to each other. Namely, we are left with, 

.\ _ p -2nm 2 K 2 G reg (x-y) 



lim d it[zi) d\ Z2) (O s {x)0- S {y) : <^iM£ 2 ) ■) = e 



lim ft ( z ,)d l f , \ 
Zi ^ z *>UiJ (2=2) 



n 2 J (fix 1 V 4 a(a;')G(x' - Zl ) - a{ Zl ) 
k 2 J d 3 x" V 4 a(x")G(x" - z 2 ) - a{z 2 ) 



In the limit r<ii and to leading divergent behavior in R, one finds 



(0 8 (x)0_ 8 (y) : (Vip) 2 (z) :> = -Ak i m 2 e- 2Km ' K ' G "'' M 



1 



(47TK)' 



-4r 



i? 4 



+ C(l/i? 3 ) + C(l/i? 2 ) 



(B9) 



b. Fusion into (V 2 <p) 2 

In a very similar set of manipulations, the fusion of two dislocation operators into the (V 2 <^) 2 operator can be 
computed. One finds the relevant quantity is given by, 

Urn V 2 Zl) V 2 Z2) (O s (x)6- 8 (y) : ^iM^) :> - e- 2 ™'"' "'^ 

Urn Vf Zi) V 2 {z2) k 2 J d 3 x' \7 4 a{x')G{x' - Zl ) - a( Zl ) 

x k 2 j d?x" V 4 a(x")G(x" - z 2 ) - a{z 2 ) 

Now, I note that V 2 z ^a( Zi ) = 2tt {5^{z l - x)9(tj — t x ) - S^ 2 \ Zi - y)B(n — t y )) which is vanishing when r <C R. 
Hence, in this case, the only terms that survive come from when the Laplacian hits G(x' — Zi ). Applying the Cauchy- 
Ricmann conditions, the fact that the Green function satisfies, V 4 G(x,i) = 4y 5^ (x)S(t) — jrid 2 G{x,t), and the fact 
that the first set of 5-functions are vanishing in the limit r <C R, an integration by parts leaves the final result, 



lim Vf Zi) V 2 {z2) (0 8 (x)0. s (y) : :) = -m A e 



Hence, one finds a familiar result. 



d t G(x - z,t x - t) - d t G(y - z,t y - r) 

(BIO) 



lim V 2 Zi) V 2 Z2) (a 8 (f)0_ 8 (y) : :> = -IGm't'e- 2 ™'"' "'^ 



e -B 2 /4T 

8ttkT 2 



(Bll) 
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c. Fusion into (Vi/p) 4 

By a similar set of manipulations, this amounts to the computation, 

hm a Uzi) ^ 2) ^ (z3) ^ 4) (6 8 (f)d_ 8 (y) : <pfr)<p&)<p&)<pfr) :) = e^™ 2 * 2 ^^) 



\im d h(zi) di 2) d Uz3) d J M 



k \ d A x'" V 4 a(x")G{x" - 23) - a(z 3 ) 



k 2 / d 3 f V 4 a(a;')G(a;' - z x ) - a{z x ) 



k / d A x" V 4 a(x")G(x" - z 2 ) - a(z 2 ) 



k \ d 6 x"" V 4 a(x")G(x" - z 4 ) - a(z 4 ) 



This is similar to the previous calculation, and one finds 

(0 8 (x)6_ s (ff) : (V#(z) :) = (4 K 4 m 2 ) 2 e~ ^^„„ (r) 



1 



2 / 2 (^) 
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